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Аннотация 

Введение. Статья посвящена проблеме исследования связанных нестационарных термоэлектроупругих полей в 
пьезокерамических конструкциях. Кратко излагаются основные подходы, связанные с построением общего 
решения при исследовании исходных несамосопряженных уравнений, описывающих рассматриваемый 
процесс. Целью работы является построение нового замкнутого решения осесимметричной задачи 
термоэлектроупругости для длинного пьезокерамического цилиндра. 

Материалы и методы. Рассматривается длинный полый цилиндр, электродированные поверхности которого 
подключены к измерительному прибору с большим входным сопротивлением. На цилиндрических 
поверхностях пластины задана температура, изменяющаяся по времени. Используется гиперболическая теория 
термоэлектроупругости Лорда-—Шульмана. Замкнутое решение построено с помощью обобщенного метода 
конечных интегральных преобразований. 

Результаты исследования. Разработанный алгоритм расчета позволяет определить напряженно- 
деформированное состояние цилиндра, его температурное и электрическое поля. Кроме того, появляется 
возможность исследовать связанность полей в пьезокерамическом цилиндре, а также проанализировать 
влияние релаксации теплового потока на рассматриваемые поля. 

Обсуждение и заключения. Использование допущений о равенстве компонент тензора температурных 
напряжений и отсутствие влияния температуры на электрическое поле позволяет сформулировать 
самосопряженную исходную систему уравнений и построить замкнутое решение. 


Ключевые слова: термоэлектроупругость, гиперболическая теория, нестационарная связанная задача, длин- 
НЫЙ пьезокерамический цилиндр, конечные интегральные преобразования. 
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Введение. В последнее время широкое распространение получили технические устройства различного 
назначения, изготовленные из пьезокерамического материала. Здесь особое место занимают приборы, работа 
которых основана на эффекте связанности упругих, электрических и температурных полей [1]. Для описания их 
работы с учетом  связанности полей в настоящее время разработаны различные теории 
термоэлектроупругости [2—4]. При этом для более качественного описания и оценки нестационарных процессов 
в конструкциях возникает необходимость построения аналитических решений. Однако математическая 
формулировка рассматриваемых задач включает систему несамосопряженных дифференциальных уравнений в 
частных производных, интегрирование которых связано с большими математическими трудностями. 

Для решения данной проблемы, как правило, исследуют уравнения в несвязанном виде [5, 6], 
анализируются бесконечно длинные тела [7-11] или рассматриваются задачи термоэлектроупругости в 
квазистатической постановке [12, 13]. 

В настоящей работе рассматривается связанная динамическая задача термоэлектроупругости для 
бесконечно длинного полого пьезокерамического цилиндра. В результате преобразования исходных расчетных 
соотношенийудается сформировать самосопряженную систему уравнений, интегрирование которой 
осуществляется методом неполного разделения переменных в виде обобщенного конечного интегрального 
преобразования [13]. 

Материалы и методы. Пусть полый, длинный, незакрепленный в радиальной плоскости, 


пьезокерамический цилиндр занимает в цилиндрической системе координат (+.,0,=) область ©: 
{а<к <Ь,0<9<2л, © << <5}. На цилиндрических поверхностях задана температура в виде следующих 


нестационарных функций (граничные условия 1-рода) — в(1.)(к=а), @>(&.) (1. =) . Внутренняя 


электродированная поверхность заземлена, а внешняя подключена к измерительному прибору с большим 
входным сопротивлением (режим электрического холостого хода). 

Математическая формулировка рассматриваемой осесимметричной задачи в безразмерной форме включает 
дифференциальные уравнения движения, электростатики и теплового баланса на основании гиперболической 
теории Лорда-Шульмана, а также краевые условия [2, 7, 15]: 


92 
ей В Ой ей уни 8 Рыб. (1) 
д’ ОГ г ог г 2 
у28 и ке СИ га, У© =0, 
дг ог г ог 
09 [2 2? 
—-| —+ О+а.УИ)=0; 
дк [2 С 20 
ви а а 
г=К,1 2; + а -- 7 9=0,9,,=®,0,„, =5,, (2) 
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1 0 т ры 
00-00 е (3) 
01 [41 


где (бр = [тра ь =, {вое} = |0", Т,,@% пер ый ы м _ в. 


33 33 


р® 


2 2 
1 71 езз _ @@зз _ &з6@зз Ь |Сз —_ Бу _ Св 
а ‚ Ч ‚ 4 =К ‚ а ‚ а 
С;зЕзз 23з37зз А\ р Л С. р С» 

И (к.,ь), 9 (..ь), 9..5) — соответственно радиальная составляющая вектора перемещений, потенциал 
электрического поля и приращение температуры в размерной форме (©^"(#.,1.) =Т(&;,1,)-Т,(.) — текущая 
температура и температура первоначального состояния тела); С.,р,е„,=. — модули упругости, плотность, 
пьезомодули и коэффициент диэлектрической проницаемости электроупругого анизотропного материала 


(", $5ЕЪ 3) ;У7з— компоненты тензора температурных напряжений (у, =Са,, у. =Сцза ); 


Л, К, а — коэффициенты теплопроводности, объемной теплоемкости и линейного температурного расширения 


материала; 8,—щ компонента тензора пирокоэффициентов; В.,—щ время релаксации; ®,— известная в 


е1 

в 9 1 
начальныи момент скорость изменения температуры; У = д НР 
гг 


В случае заземления внутренней поверхности пьезокерамического элемента, электрическое напряжение 


У(.) определяется потенциалом на его внешней поверхности: 


У(1.) =9(1,1.). (4) 
При построении общего решения на первом этапе исследования в результате интегрирования уравнения 
электростатики определяется радиальная компонента вектора напряженности электрического поля: 


ми 0 р 


т г. (5) 


где р, — постоянная интегрирования. 
Подстановка (5) в (1)-(3) позволяет сформулировать новую задачу относительно функций И (",1),© (1,1) ; 


При этом условие отсутствия радиальной составляющей вектора индукции электрического поля на внешней 
цилиндрической поверхности элемента (последнее равенство (2) выполняется в случае 2 =0, а условие 


заземления внутренней поверхности ( ф,_, =0) удовлетворяется в результате интегрирования (5). 


На следующем этапе решения неоднородные краевые условия (2) приводятся к виду, позволяющему в 
дальнейшем использовать процедуру неполного разделения переменных методом конечных интегральных 


преобразований. Для этого вводятся новые функции и(г,!), №(г,!) связанные сИ (г,!),©(г,!): 
О (т.г) = Н, (г.р) -+и(т,г),©(г,г)= Н, (г,1)+М (т, г), (6) 
где Н, (и,г) = Л (г) А(1,1) + (К) А(К,1)+ № (г) <, (1) + Л, (г) (1), 
Н, (гг) = № (г) (1) + , (г) < (1), 
12)... 1,(г) — дважды дифференцируемые функции, А(т.г)_„, = (1+а, -а; а, И (т.г). 
Подстановка (6) в расчетные соотношения (1)-(3) относительно функций И (",!),©(т,г) при выполнении условий: 
(1+а, )УН, +(а, -ЮН, =А(,#) Нк =@, Ну =4, (1) 


г=ЮЛ ? 
позволяет получить новую краевую задачу относительно функции и (, г) р М№(", 1) : 


ди и ОМ № 1 2би м 


У р +В +В Е 8 
бе ‘о бд Ак (1+а,) 0’ ' - 
ОМ 1 5? 
У -[бьр-2 аумчаун) = 5 
г= Ю,1 Уи=0 , М=0: (9) 


=би=-Н, (.0),М=-Н, (,0), (10) 
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би _ СН, (г,1) М 6 @Н, (г,г)_ 
61 в м ^ $ 
О а. 
где Е, =-У = [Сев ен, +а,УН,), 
2 
Е ей е +Ь, р +Ь, -ь - — 
дг г дг г (1+а,) м 
ь _ а, +а.а, ь _ @&-1 в -ь+ 26 
1+а, Г 1+а, . 1+а, 


Здесь следует отметить, что А(х Г) является функцией перемещений цилиндрических поверхностей цилиндра. 


Первоначально А(т, 1) приравнивается к нулю с последующим ее определением и уточнением Н,,Ё,Ё. 
Дальнейшие преобразования расчетных соотношений (8)-(10) связаны с использованием следующих 
допущений: В, =1,Б, =0 и введения термоупругого потенциала 
№=УВ. (11) 
Условие В =1 можно принять без большой погрешности, поскольку для пьезокерамических материалов 
Ь =0,94--0,98, а зависимость В. =О выполняется в случае равенства компонент тензора температурных 


напряжений (у, =У., ) и отсутствие влияния температуры на электрическое поле ( з. =0). 


В результате формируется следующая задача относительно и 1 В В (. Г) : 


д д 1 2и 
— Уи-+Ь УВ =Ё, 12 
ОГ > бк (1+а, ) ОР - ие 
д 12) 27 
був в. вчнии)= К; 
г=А,1 Уи =УВ =0 ; (13) 
=Ои=-Н, (,0) ‚УВ=-Н, (*,0}, (14) 
би _ СН, (г,1) 72) ув-@,—2? (^,1) | 
Ге Ге Ге Ге 
ОН 12) 5? 
где Е, =- = +2 ое, +а,Н,),УН, = № (г)@ (1) № (г) ®, (1). 


Начально-краевую задачу (12)—(14) решаем, используя структурный алгоритм обобщенного конечного 
интегрального преобразования (КИП) [14]. При этом для данной задачи удается использовать 


однокомпонентное неизвестное ядро преобразований К (^, г) } 


{С , С, (1) = ии) вк} рьтак, (15) 


{и(т,г),В(к.г)} = >16 (д.0, са, вук (лк, (16) 


: 
К, | — [ка, к)? ; 
К 


где Я, — собственные значения, образующие счетное множество. 


В результате использования алгоритма КИП [14] получаем задачи относительно ядра преобразований К (1, з г) : 


"КА К (Я, 
4 а ( -?-( ых Ка,=0. (17) 
аг г @ г 
Г ЕЛ УК(А,—)=0, (18) 
и трансформант С (4,,1),С.(4,,Ю: 
Я: 1 4С, 
С, -Е - “ =Ён › (19) 


(1+а,) “ (1+а,) а 
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2 А 4? | 
Я, С, - а (а.б,+а,С,)=Ен , 


АС", ас, 
=” бы =би, —“=Си; (20) 


=0С,=С 
ИИ 


10 › Я 


где{Ен, Ен} = [{Е ‚Е } К(Аьгтаг, {Си 6} =-|{Н, (т,0),Н, (›,0)} КС „гта, 


4 в Г ©Н, (г.г) | ЯН, (+) 
{Си» бью} «о аа 5 К(А,,гугаг. 


Общее решение задачи (17), (18) имеет вид: 
К(А,") =\ (4). Л, (Яг)- Ло (4) (4). 20 
Здесь собственные значения 7, определяются с помощью следующего трансцендентного уравнения: 
% (а) (Ав) - (4) (Ав) =0. 
Система дифференциальных уравнений (19) приводится к следующему разрешающему уравнению 4-го 


порядка относительно С,(,,№): 


И а, м 
[2 +, аг +Ь, аг +В Р +, |6; =Рн, (22) 


ог 


1 Я: 
ира : | ! ыы 
В а’ ар В а, а. В 


Поскольку характеристическое уравнение, соответствующее (22), 
КНЬ К +Ь, К? +В К-Ь,, =0, 


(1+ ь 
где Е, Ав, +(1+а,) к,|- 2 ый С к», 
1 


будет действительное, то оно из условия осциллирующего решения для С, (1,1) имеет два действительных 


корня (А, К., ) и два комплексно сопряженных корня 
(К; ЕС, К, =9-№). 
В этом случае общее решение уравнения (22) имеет вид: 


С, (4,1) = Бы ехр(ки)+ Бу, ехр(к.1)+Р., ехр(К;/) + Бу, ехр(К. 1) (23) 


+, (г) [ехр [^, (1 -т)] —ехр[ А, (:— г)] Ат [в (г)ехр[ ©(1 -- г) | (6 зт(а/—@т)-Ь, с0$(@# + @т)} Ат, 


где, = [(, К ) (®, -в)} +07] в = (+5, а ы к. Бо; = КиК (К, +) 0+0? —@°, 


—К 
2 


К, 41 


В: ==): @= 


Функция С,(Й,,Г) определяется из первого уравнения системы (19). Подстановка полученных выражений 
для трансформант в граничные условия (20) позволяет определить постоянные интегрирования Д,,...О.,. 
Подстановка С1(Я,,1) ‚С. (Я,,В)в (16), (11), (6) позволяет получить окончательные выражения для функций 


О (+, 1) ‚@(х»,2): 
(и) = н, (+ а,Ок (АКГ. 64) 


@(,.:) = Н, (.1)+У ба ,бУк (ак 


11 
На заключительном этапе исследования функции Н, (+,г),Н, (1,1) определяются при решении следующих 


дифференциальных уравнений: 
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ие 5; (25) 
дг Г ог 


что позволяет существенно упростить правые части (Е, Е) расчетных соотношений (8). 


Подстановка выражений для Н,,Ньв (25) позволяет сформировать системы уравнений относительно 
функций 1 (1)... (г) ‚ которые определяются при удовлетворении условий (7). 

Потенциал электрического поля пьезокерамического цилиндра определяется в результате интегрирования 
равенства (5) и удовлетворение предпоследнего граничного условия (2): 


ОН, (г,Е Н! (г 
#- | ц м г. ) на.н, бы (26) 
+в „ов“ +4, У. С, (А.В, (4) [к +10% 4) 
1 1=1 
Н Н о — 
пая ии) вождю 
Г Г = 
Р _ 2К(А, К(4, 
+, »`6, (4. ВВ (ЯК, | й ‚ В, (4-1 и и и 2. 
1=1 |\=А 


В, (А) = [УК (4. 


Полученные расчетные соотношений (24), (26) удовлетворяют дифференциальные уравнения (1) и краевые 
условия (2), (3), т.е. являются замкнутым решением рассматриваемой задачи. 

Результаты исследования. В качестве примера рассматривается радиально поляризованный 
пьезокерамический цилиндр (Р=0,02м, К =0,8) состава РИТ-4, имеющего следующие физические 


характеристики [10]: р = 7500 кг/м °, 


{С»С„Сз} = 113,9; 11,5; 7,43} ж10® нм? {е,»е,} = [-5,2; 15,1} Клум”, 
{уз} = 4,6; 3,9} <10°НИм? °С), =, = 5,6210 фм, а, =2х10 *Кл/м? °С), 
К =3х10° Дж/(м * °С), Л =1,6 Вт/(м °С), В. =5х10с. 


На внутренней поверхности ( 7, =а ) пьезокерамического цилиндра действует температурная нагрузка: 


а) «57 ес. кн.) ‚а (1,)=0, 


тах 


где Н (Г) — единичная функция Хэвисайда ( Н (7) =1 при Г>0, Н(Р)=0 при #<0), Ти =Ти, 1, Тик 
— максимальное значение внешнего температурного воздействия и соответствующее ему время в размерной 
форме (Т", = 373 К (100 °С), 7, =293 К (20°С), #'„ =1 0). 

На рис. 1 представлены графики изменения функций ©” (г,г), ( (т,г), $(г,г) по радиальной координате г в 


различные моменты времени 1. Цифрами 1-3 соответственно обозначены результаты, полученные при 


Л * 
151 


тпах тах (ах — Кр? ах ). 


следующих значениях времени: #=1 41 


тах ? 


Анализ результатов расчета позволяет сделать следующие выводы: 

— достаточно большая величина коэффициента линейного температурного расширения @, пьезокерамического 
материала приводит к быстрому прогреву цилиндра; 

— радиальные перемещения на внутренней цилиндрической поверхности ( 7 = К) на первом этапе исследования 


(1 ЕР ) принимают наибольшие значения с последующим уменьшением в течение времени. Обратная картина 


тах 


наблюдается относительно перемещений при г=1; 
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9° (+,:), °С 


80 


60 


40 


Рис. 1. Графики изменения ©" (+, 1 й и (+, УФ ("1 ) по радиальной координате г в различные моменты времени # 


(2-4 щ›3—154 15 ): а— ©" (и,Р) к; 6—0 (и) -+г;в— ф(и,4)+к 


Степень связанности термоэлектроупругих полей наиболее удобно проанализировать с ПОМОЩЬЮ 


Я? аз й 
коэффициента Б,, = т 1+а, +— | равенства (21). Здесь а, определяет связанность электроупругих полей, а 
Я 


а. 
— влияние скорости изменения объема тела на его температурное поле. 
а 

7 


На рис. 2 представлен график изменения перемещений И (1, 1) по времени { с учетом (сплошная линия) и 


без учета (пунктирная линия) наведенного электрического поля. 
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Следует отметить, что предварительная поляризация пьезокерамики приводит к образованию более 
«жесткого» материала (а, =0,353) и соответственно к уменьшению перемещений при деформировании 
цилиндра. 

Связанностью температурного и электроупругого полей в пьезокерамическом цилиндре можно пренебречь 
за счет малого значения 


48 18х10 < 1. 
а 


(г, Г) х10’ 
10 


2х10° 3х10° 


Рис. 2. График изменения И (1, 1 ) по времени 1 


(сплошная линия — с учетом поляризации, пунктирная — без учета электрического поля) 


На рис. 3 показаны графики изменения электрического напряжения У (1) по времени с учетом (сплошная 


линия) и без учета (пунктирная линия, р, =0) релаксации теплового потока. 


У х!0’ 


0 100 200 300 


Рис. 3. График изменения У (1 ) по времени # (1. „ =0,001 с,Ё„ =196) 


Результаты расчета показывают, что для рассматриваемой задачи уточненную гиперболическую теорию 
Лорда-Шульмана необходимо использовать при большой скорости изменения температурной нагрузки 


ав (1.) 
Ги 


(г ‹ 20,001 с, >5,56х10° к ), а при меныших скоростях — классическую теорию 


* |" =0 
термоэлектроупругости (Д., =0). 

Обсуждение и заключения. Построенное новое замкнутое решение связанной динамической задачи при 
удовлетворении граничных условий теплопроводности 1|-го рода позволяет определить все компоненты 
термоэлектроупругих полей в длинном пьезокерамическом цилиндре. Преимущество представленного 
алгоритма расчета заключается в том, что в отличие от несвязанной постановки задачи, отпадает 
необходимость аппроксимации функции температуры при исследовании уравнения движения. При этом 
действительно влиянием скорости изменения объема пьезокерамического тела на его температурное поле 
можно пренебречь. 
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